Collision-Dependent Atom Tunneling Rate in Bose-Einstein Condensates 
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We show that the interaction (cross-colhsion) between atoms trapped in distinct sites of a double- 
weU potential can significantly increase the atom tunneling rate for special trap configurations 
leading to an effective linear Rabi regime of population oscillation between the trap wells. The 
inclusion of cross-coUisional effects significantly extends the validity of the two-mode model approach 
allowing it to be alternatively employed to explain the recently observed increase of tunneling rates 
due to nonlinear interactions. 
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Atom optics and the physics of ultracold matter 
waves have witnessed rapid theoretical and experimen- 
tal progress since the achievement of atomic vapor Bose- 
Einstein condensates (BECs) The interest in such a 
system is quite wide ranged as it has opened new tech- 
nological frontiers 0, 1^, and has renewed the investi- 
gation on many-body physics fundamental issues once 
BECs trapped in optical lattices offers a powerful toolbox 
for implementation of condensed matter model systems 
with highly controllable parameters H, |^ . 

The simplest system to exhibit non-trivial many-body 
phenomena is the two-mode BEC trapped in a double 
well potential. For this system it was previously pre- 
dicted 0i a plethora of dynamical regimes, including 
Josephson oscillation (2| of atoms between distinct trap 
sites, and macroscopic self-trapping of atomic popula- 
tions, where coherent tunneling is suppressed when the 
number of trapped atoms exceeds a critical value. The 
first realization of an atomic single BEC Josephson junc- 
tion has indeed experimentally demonstrated ^3 those 
predictions Moreover an increase of one order of 

magnitude of the atomic tunneling rate (from 2 to 25 Hz) 
due to nonlinear atomic interactions was then observed, 
which cannot be explained inside those previous models. 
In those models 0, B a two-mode approximation (TMA) 
is invariably applied, in which each trap well is popu- 
lated by one single localized condensate mode. However 
due to a localization assumption the BEC dynamics is 
considered by regarding only self-collisions (on-site inter- 
action) between atoms of individual condensate modes. 
Notwithstanding, distinct trap configurations can be as 
such that cross-collisions (off-site interaction) of atoms 
in distinct BEC modes cannot be neglected and the ob- 
served coherent phenomena are affected by cross-collision 
induced effects. This situation has been previously con- 
sidered either by directly integrating a non-polynomial 
Schrodinger equation [ifll, or by extending the TMA 
through Gross-Pitaeviskii equations [l^. 

In this paper we investigate the cross-collisional ef- 
fects over the quantum dynamics of a BEC trapped 
in a double-well potential inside the TMA. By appro- 
priate manipulation of the trapping potential the sys- 



tem shows a very rich quantum dynamics related to the 
cross-collisional rates. In absence of cross-collisions, self- 
trapping is observed in agreement with previous theo- 
retical predictions 0, 0|. Remarkably, the inclusion of 
cross-collisional effects can inhibit self-trapping, result- 
ing in an effective coherent oscillation (effective Rabi 
regime) of atoms between the two wells with an effec- 
tive Rabi frequency, dependent on the cross-collision rate 
and the number of atoms in the sample. Moreover cross- 
collisions increase the range of validity of the TMA (when 
the many-body interactions produce only slight changes 
on the ground-state properties of the individual poten- 
tial sites). The new limit can be as such that the TMA 
is indeed valid for the number of atoms present in the 
experiment of Ref. [Tflj , an d thus can be employed for 
alternative explanation Il2l | of the observed increase 
of tunelling rate through off-site interactions. 

We consider an atomic BEC trapped in a symmetric 
double-well potential V{r) with minima at ri, and r2, 
such that V^(ri,2) = 0, the many-body Hamiltonian is 
given by 

fc2 
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where m is the atomic mass, Uq = Anh^a/m measures 
the strength of the two-body interaction, a is the s-wave 
scattering length, tp{r,t) and ■0^(r,i) are the Heisen- 
berg picture field operators. We must emphasize that 
throughout our calculations we assume that the bar- 
rier height separating the two wells is sufficiently larger 
than the condensate chemical potential, Vq > iic, so 
that a TMA can be considered to describe the system 
|l3l| . The derivation of the two- mode Hamiltonian in 
a double-well trap follows the standard procedures of 
Ref. - The lower energy eigenstates of the global 
double- well are approximated as the symmetric and anti- 
symmetric combinations u±{r) w -^[ui{r) ±it2(r)], and 
the field operators are expanded in terms of the local 
modes Uj{r), j = 1,2 and the Heisenberg picture an- 
nihilation operator as aj{t) — J d^ru*{r)il;{r,t) so that 
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[aj,a|,] ~ Sjk- The potential expanded around each 
minimum is V{r) = V'-^^r — rj) + ... j=l,2, where 
\/(2)^j. _ j.j'j jg harmonic approximation to the po- 
tential in the vicinity of each minimum. The normalized 
single-particle ground-state Uj{r) of the local potential 
V(^)(r-rj), with energy Eq , defines the local mode solu- 
tions of the individual wells. The tunneling frequency H. 
between the two minima is then given by the energy level 
splitting of these two lowest states, = 2R/h, where 
R — J d'^rul{r)[V{r) — V'--^\r — ri)]u2(r). By noticing 
that the total number of atoms is a conserved quantity, 
after some algebraic manipulation, the following Hamil- 
tonian is then obtained in the interaction picture 

H = n[2A{N - 1) + n]{ala2 + alai) + hr]{a\a2 + alai)^ 
+n{K - Tl)[{a\f{al) + {a\f{al)] + nr^N{N - 2) (2) 



where 



~ 2hV^ 



If 



is the self-collision rate and rj = 



(II )/ d^ru*UjU*UQj, and A = i^) J d^vu*UjU*Ui, are 
cross-coUisional rates. Here V^^j = / d^rlujl"* is the j- 
mode effective volume. It is immediate that for fin — > 0, 
the Hamiltonian |2l reduces to that previously investi- 
gated in the study of tunneling in condensates 0]. Also 
for ?i?7 — !■ and hK — > 0, the Hamiltonian ^ reduces 
itself to that discussed in 0, Q , accounting for tunnel- 
ing oscillations as well as for population self-trapping. 
As 77 ^ K in Hamiltonian ^ a new dynamical regime 
of stable long-time tunneling (hereafter called effective 
Rabi regime) can be attained. 

Introducing the Schwinger angular momentum opera- 
tors Jj, — (d\a2 — a|di)/2, Jj, = i{a\di — a|a2)/2, = 
{a\d2 + d\ai) /2, the Hamiltonian (0) can be rewritten as 



(3) 



where we have neglected constant energy shifts depen- 
dent on A'' and have defined a new effective tunneling 
rate W = 2[U + 2A{N — 1)], which is dependent on the 
number of atoms in the atomic sample and the cross- 
coUisional rate A. In fact A is very small compared to SI, 
but for a sufficiently large number of atoms the additional 
tunneling term may lead to observable effects whenever 
2A{N - 1) >n. Furthermore the third term of Eq. © 
shows that the cross-collisional rate 77 competes with the 
self-collision rate k leading to an effective on-site colli- 
sion rate k' = k — r]. Since < 77 < k, k' on Eq. Q 
can be disregarded in a trap configuration where i] ^ k 
with f2' ^ k' . To infer the validity and implications of 
such a regime we consider a specific trapping potential 

[3 of the form V{r) ^ b {x^ - j^f + ^muj^iv^ + z^), 
where the inter-well coupling occurs along x, and LUt is the 
trap frequency in the y-z plane. This potential has fixed 
points at = The position uncertainty for a har- 
monic oscillator in the ground state is xo = \ / — , with 
ojQ = ^ Ad/m. For a suitable choice of the barrier height 



only two energy eigenstates lie beneath the barrier. As- 
suming LOt = OJQ the local mode on each well is then given 
by ". (^) = il^Y^ exp (z[(fiZ^^±^) . The collision 
rates may then be evaluated to give k = i 
77 — Kexp(— (70/2x0), and A = k exp(— 3gQ/4xQ), while 
17 = 2^£exp(-goV2x§). 

Firstly we shall investigate the validity of the localiza- 
tion of the TMA considered in the present description. 
It is valid only if the many-body interactions produce 
small modifications on the ground-state of the individual 
potential wells 15]. In the absence of cross-collisions 



the TMA is fimitcd to ?icJo = 



> NhK = N 



n\UQ\ 



and thus the BEC is limited to a few 100 atoms. The 
inclusion of cross-collision necessarily increases the effec- 
tive mode volume, which in terms of the effective on-site 
collision rate n' is given by V^^^ = V'e///(1 — e"''"/^^"). 



For the trapping potential considered, Veff ^ Stt'^/^; 
and the TMA is valid whenever 

2n^/^x^ 



0' 



N <^ 



\a\{l - e-io/^K) 



(4) 



For suitable values of the trap parameters, it is thus pos- 
sible to have very large numbers at the right of the in- 
equality above. Hence the considered TMA, with cross- 
collision included, can be appropriately employed for the 
description of the observed phenomena in present exper- 
imental setups 10] for numbers of atoms N > 10'^. 

Typically 77, A ^ SI, but cross-collisional effects can 
only be neglected if 2k{N ~ < 1, i. e. if 

{N - 1) °/2°2 e"^"^'^'^" < 1- By considering the exper- 
imental values from Ref. 10] for = 1300, — 5.3 nm 
for ^^Rb atoms, and go = 2.2/im, we obtain that only for 
xq ^ 0.82 go is that cross-collisional effects can be ne- 
glected. This bound value however depends strongly on 
the number of atoms in the sample and decreases rapidly 
as A^ is increased. For example, if the number of atoms 
is increased in one order of magnitude to TV = 13000 
then the cross-collisional effects can be neglected only if 
xq <^ 0.366 go. Furthermore by equating SI' to the ob- 
served frequency of oscillation of 25 Hz in the experiment 
performed by Albiez et al. ^3 and Q. to the expected fre- 
quency of tunneling of 2 Hz we find with the above pa- 
rameters that only for xq « 6.56 go is that the tunneling 
frequency would significantly be increased due to non- 
linear off-site interactions. Of course, inside the applied 
TMA this situation corresponds to that of two strongly 
overlapped modes, which would compromise the validity 
of the localization assumption and the considered gaus- 
sian eigenfunctions. However the relation between xq and 
go is not linear and depends strongly on TV and on go as 
well. For a fixed go the xq to go ratio decreases rapidly 
as A^ is increased. Again, for A^ = 13000 the substantial 
increase in the tunneling frequency would be observed for 
Xo ~ 0.89 go. We remark that the above analysis could 
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be similarly applied by considering the local modes as ap- 
proximately given by the excited states of the harmonic 
oscillator. In such a case the localization conditions in 
relation to the distance between minima of the trapping 
potential would be significantly relaxed |18l] . 
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FIG. 1: Numerical results for the eigenvalues of Eq. 3 in the 
eigenbasis of Jx- The energy spectra is shown for different 
regimes of the ti/k ratio. Each spectrum has been displaced 
by a constant factor. 

We now turn to the description of the dynamical 
regimes attained by the Hamiltonian (PJ. In Fig. 1 
we plot the numerically calculated eigenvalues spectra 
for a condensate with 1000 atoms and fixed JI/k = 50 
ratio. Each curve in the plot represents a specific rj/n 
ratio as depicted in the figure. For ti/k — 1/1000 it is 
immediately verified that the eigenvalues correspond to 
those of the model without cross-collision discussed in 
Ref. 0- The low-excited eigenstates are close to the 
eigenstates of the Jz operator given by |j, — j)z- Highly 
excited eigenstates are closer to the eigenstates of the 
operator J^, given approximately by |j, O)^;. This last 
state suggests an ideal situation for occurrence of popu- 
lation self-trapping for large N . The infiexion point for 
intermediate eigenvalues represents the frontier for the 
occurrence of self-trapping. As the ratio n/rj is increased 
the eigenvalues behavior changes dramatically due to the 
influence of the term proportional to An outstand- 
ing feature is that the inflexion point is shifted to higher 
eigenvalues for k/j] = 1/100 — 1/10 and completely disap- 
pears for k/t] > 1/5. For k/t] w 1 the term proportional 
to is negligible in comparison to the other terms and 
the Hamiltonian eigenstates get as closer to the Jz eigen- 
states as k/?7 1. 

The equations of motion for the angular momentum 
operators are given by 

Jx = -hfl' Jy ~ Ahr][Jy, Jz]+ (5) 
Jy ^ hn'X -2h{K-3T])[Jz,.U+ (6) 
Jz = 2h{K-r])[Jy,X]+ (7) 
One can immediately see that if the cross-collision term 



is strong enough so that k — rj <S^n' then 

Jz « (8) 

and thus Jz being a constant of motion the set of Eqs. 
El and IS)) is solved for J^ as 

(t) = 4(0) cos n" t + Jy (0) sin n" t, (9) 

where rj" = W + 8r]Jz = [n + 2A{N - 1)] + SrjJz is the 
new Rabi frequency, which now depends on both the 
number of atoms and the initial condition for Jz through 
the nonlinear cross-collision rates A and 77, respectively. 
Since rj <^ ^l, Jz value only gives a small shift in f2 ' . 
The effective Rabi regime of 77 ^ k is justified whenever 

In Fig. 2 the time evolution of the mean value (Jx) 
as given by the general system of Eqs. (5-7) for distinct 
tj/k ratio. We suppose an initial state localized in one 
well as \j,—j)x- Similarly to Ref. due to intrinsic 
quantum fluctuations in the initial condition there are 
some oscillations of the quantum mean decay |l9l| . In 
Fig. 2(a) the revival of the oscillation that occurs at 
later times is again due to the discrete spectrum of the 
many-body Hamiltonian [20I l2l|. In Fig. 2(b) the self- 
trapping of the atomic population is shown. This result 
should be compared to that of Fig. 2 from R.ef.[lol|. More 
interesting though are the new features in the quantum 
dynamics due entirely to the presence of cross-collisions 
appearing in Fig. 2(c). Notice that even when cross- 
collisions are relatively small {n/rj ^ 10) coherent tunnel- 
ing dominates, suppressing the initial self-trapping. As 
the cross-collisions become stronger the collapse and re- 
vival dynamics changes. For large intervals of time there 
is a total collapse of the mean value, which lasts for large 
periods of time in such a way that this period decreases 
as 77 gets closer to k approaching the limit where it tends 
to zero when the effective Rabi regime takes place and 
there is no collapse and revival. Notice that the decreas- 
ing amplitude modulation in the collapse and revival is 
due to nonlinearities introduced in Eqs. (5-7) by cross- 
collisions. As 77 — !• K however, the amplitude of coherent 
oscillation of the population difference goes to a constant 
value with oscillating frequency f2 . 

To push the discussion presented here further, it is in- 
teresting to consider the impact of the cross-coUisional 
terms in related many wells systems 0. The dynamics 
is richer and additional effects may occur as consequence 
of the overlap of neighboring condensed modes. The re- 
sults found for the double- well case can be extended and 
set into the more general framework of hopping in many 
wells systems. It is possible to render a general discus- 
sion of the cross-coUisional terms for the derivation of 
a discrete extended Bose-Hubbard model 0| and of the 
Wanier-Stark effect based solely on ground state trap 
solutions. For instance, a 77-mode ansatz of the natural 
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extension of the Hamiltonian in Eq.(2) results in a n- 
mode extended Bose-Hubbard Hamiltonian 

H = '^{ni - 'nLi+i)niini - 1) 

+ J2m,i+i{al+iai + H.cf (10) 

I 

with all parameters related to the those for the double- 
well system. The dynamical regimes given by (|10|) have 
in part a similar aspect to the ones presented for the 
double-well potential but can be extended since all the 
on-site and off-site parameters can be independently var- 
ied. Those issues are being further investigated in the 
context of phase transitions and shall be addressed else- 
where Ei. 



ity of such a two-mode approximation increases within 
the assumption of strong cross-collisions in such a way 
that it may be in accordance with experimental data [lOj 
observed for the increase of tunneling frequency due to 
nonlinear interaction. As a last result the collapse and 
revival dynamics is changed resulting in large periods of 
total collapse. As r] ^ k, however the collapse and re- 
vival frequency increases attaining the regime of Rabi 
population oscillation for rj ^ k. It is remarkable that 
similar results to tho se p reviously given in exact numer- 
ical calculations [10, llJ could be derived inside the 
TMA just by adding cross-coUisional effects. We hope 
that the above considerations bring some contribution to 
experimental investigation. 
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FIG. 2: Numerical results for the mean value of the popula- 
tion imbalance {X) for N=100, kN = 2.0, (a) r? = k/10, (b) 
77 — k/100 and {c)ri = k/2. For convenience we normalize the 
time in units of the single particle tunneling period so that 
in these units Q — 1. Time evolution found by integrating 
Schrodinger equation in the basis for the Hamiltonian in 
Eq. 3. 

In summary we employed the TMA to model an atomic 
two-mode BEC trapped in a double well potential con- 
sidering cross-collisions between the condensed modes. 
Cross-collisions strongly inhibit the self-trapping phe- 
nomenon even for small values of the cross-collision rate 
ry. For a given trap geometry the eigenvalues of the many- 
body Hamiltonian determine a new dynamics resulting in 
an effective coherent oscillation (effective Rabi regime) 
of the population in both wells. In the limit r; — > k for 
K — 77 <^ 17 an effective Rabi regime takes place with an 
effective Rabi frequency, f2" , which is explicitly depen- 
dent on the total number of atoms. The limit of valid- 
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